Unconventional motional narrowing in the optical spectrum of a semiconductor 

quantum dot 
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Motional narrowing refers to the striking phenomenon where the resonance line of a system cou- 
pled to a reservoir becomes narrower when increasing the reservoir fluctuation. A textbook example 
is found in nuclear magnetic resonance, where the fluctuating local magnetic fields created by ran- 
domly oriented nuclear spins are averaged when the motion of the nuclei is thermally activated. 
The existence of a motional narrowing effect in the optical response of semiconductor quantum dots 
remains so far unexplored. This effect may be important in this instance since the decoherence 
dynamics is a central issue for the implementation of quantum information processing based on 
quantum dots. Here we report on the experimental evidence of motional narrowing in the opti- 
cal spectrum of a semiconductor quantum dot broadened by the spectral diffusion phenomenon. 
Surprisingly, motional narrowing is achieved when decreasing incident power or temperature, in 
contrast with the standard phenomenology observed for nuclear magnetic resonance. 

PACS numbers: 78.67.Hc, 78.55.Cr, 05.40.-a 
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In the seminal work on motional narrowing by Bloem- 
bergen et al., relaxation effects in nuclear magnetic reso- 
nance were beautifully explained by taking into account 
the influence of the thermal motion of the magnetic nu- 
clei upon the spin-spin interaction [l| . The general treat- 
ment of relaxation processes for a system interacting with 
a reservoir was later formulated by Kubo in a stochastic 
theory that assumes random perturbations of the sys- 
tem by a fluctuating environment 0]. Depending on 
the relative magnitude of the spectral modulation ampli- 
tude and the inverse of the modulation correlation time, 
the relaxation dynamics is either in the slow modulation 
limit, where the optical line-shape reflects directly the 
statistical distribution of the different system energies, 
or in the fast modulation limit where the fluctuation is 
smoothed out and the line-shape is motionally narrowed 
into a Lorentzian profile. The relevance of motional nar- 
rowing for the description of relaxation phenomena has 
spread throughout many different fields, such as spin re- 
laxation in semiconductors vibrational dephasing in 
molecular physics Q, or phaTe noise in optical pumping 

a. 

The optical spectrum of a material system with local- 
ized, zero-dimensional electronic states provides a generic 
example of the influence of a fluctuating environment on 
the coherence relaxation dynamics. In that case, the per- 
turbing interactions induce a stochastic shift over time of 
the optical spectrum, resulting in the so-called spectral 
diffusion effect, which was observed for rare-earth ions , 
molecules Q , or semiconductor quantum dots 0, 13 ■ ^ n 
this latter system, impurities, defects or localized charges 
in the vicinity of a quantum dot induce micro-electric 
fields that shift the quantum dot emission line through 
the quantum confined Stark effect. The fluctuation of the 
quantum dot environment thus randomize the emission 
energy over a spectral range S on a characteristic time 



scale r c . Spectral diffusion under light illumination was 
reported with jitters of the quantum dot emission energy 
from hundreds fieV to few meV on time scales ranging 
from milliseconds to minutes 0,111 El- However, to the 
best of our knowledge, no evidence for motional narrow- 
ing (Er c <C?i) was pointed out for the spectral diffusion 
phenomenon in quantum dots. 

We present here the experimental evidence for mo- 
tional narrowing in the optical spectrum of a semicon- 
ductor quantum dot. High-resolution Fourier-transform 
spectroscopy performed on the photoluminescence signal 
of a single InAs/GaAs quantum dot allows the deter- 
mination of both width and shape of the emission line. 
A crossover from Lorentzian to Gaussian profiles of the 
zero-phonon line is observed while its lincwidth increases. 
We obtain a quantitative agreement with Monte Carlo 
simulations where the fluctuating environment is mod- 
elled by the random capture and escape of carriers in 
traps located in the quantum dot vicinity. We show that 
motional narrowing is achieved when decreasing the tem- 
perature or the incident power because of the asymmetry 
of the capture and escape mechanisms. Our study pro- 
vides a counterintuitive example of motional narrowing 
compared to the standard phenomenology described in 
nuclear magnetic resonance. 

We study self-assembled InAs/GaAs quantum dots 
grown by molecular beam epitaxy in the Stranski- 
Krastanow mode, and micro-photoluminescence mea- 
surements under non-resonant excitation are performed 
in the far field using the experimental setup described in 
Ref. I n order to accurately characterize the quan- 

tum dot emission spectrum, the photoluminescence sig- 
nal arising from a single quantum dot is analyzed by 
high- resolution Fourier-transform spectroscopy |ll| . The 
Fourier-transform technique is implemented in the de- 
tection part of the setup where the photoluminescence 
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signal passes through a Michelson interferometer placed 
in front of a 32 cm grating spectrometer. The sig- 
nal is detected by a low noise Si-based photon count- 
ing module. With a translation stage, we vary the time 
t for propagation in one arm of the interferometer and 
record interferograms of the photoluminescence emission 
I(t)=Io(l+C(t) cos(E t/h)), where Jo is the average pho- 
toluminescence signal intensity, Eq the central detection 
energy, and C(t) the interference contrast. In the two 
cases of Lorentzian or Gaussian line-shapes, a Fourier- 
transform analysis shows that the interference contrast 
decay is exponential or Gaussian, respectively. As we 
will see below, the implementation of this technique in 
single quantum dot spectroscopy allows an accurate de- 
termination of both width and shape of the emission line. 
As far as the linewidth is concerned, sub-/zeV measure- 
ments can be performed by using this interferometric 
correlation technique ^l|. For the determination of the 
line-shape, our technique allows a high-resolution sam- 
pling of the Fourier-transform of the spectrum on typi- 
cally thousands of points. This value is by two orders of 
magnitude larger than the average number of illuminated 
pixels in a charge-coupled device in the case of a multi- 
channel detection in the spectral domain. The drawbacks 
of our Fourier-transform technique are the long acquisi- 
tion times of roughly one hour per interferogram, and 
the required extremely high mechanical stability of the 
system during the interferogram acquisition. 
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In Fig. ^ we display the measured (squares) interfer- 
ence contrast C(t) for the emission spectrum of a sin- 
gle InAs/GaAs quantum dot at 10K, on semi-logarithmic 
plots, for three different incident powers: 0.18 (a), 0.72 
(b), and 2.88 kW.cm" 2 (c). We first observe that the 
coherence relaxation dynamics becomes faster when in- 
creasing the incident power. Moreover, we notice a grad- 
ual modification in the shape of C(t). At low power 
(Fig. da)), the interference contrast decay is quasi- 
exponential with a time constant of 29 ps, thus corre- 
sponding to a quasi-Lorentzian profile with a full width 
at half maximum (FWHM) of 45 /zeV. At higher power 
(Fig. Gib)), C{t) has a Gaussian decay at short times 
(i< 15 ps), and remains exponential at longer times. The 
line has an intermediate profile, and here its FWHM is 
105 /ieV. At the highest incident power (Fig. Ofc)), the 
interference contrast decay is predominantly Gaussian, 
thus corresponding to a quasi-Gaussian profile, with a 
FWHM of 155 ^eV. 

This crossover from a Lorentzian to a Gaussian line- 
shape when increasing the incident power reveals an orig- 
inal relaxation dynamics in semiconductor physics. In 
fact, in bulk semiconductors or in quantum wells, the 
standard phenomenology is the exact opposite. An in- 
creasing number of photo-created carriers leads to the 
activation of carrier-carrier Coulomb correlations, and 
induces the inverse transition from an inhomogeneously 
broadened Gaussian line to a homogeneous Lorentzian 
one |12| . In the following, we provide a quantitative anal- 
ysis of our measurements in the framework of the Kubo 
theory. 
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FIG. 1: Interferogram contrast C(t) of the photoluminescence 
signal of a single InAs/GaAs quantum dot at 10K, on semi- 
logarithmic plots, for three different incident powers: 0.18 
(a), 0.72 (b), and 2.88 kW.cm -2 (c). Data (squares), system 
response function (dotted line, in red), theoretical fits (solid 
line, in green) obtained by the convolution of the system re- 
sponse function with Eq. Q, with Er c /ft~ 0.6 (a), 1.05 (b), 
and 1.35 (c), are plotted as a function of the delay t. 



exponential decoherence (C(i)~exp(— I] 2 T c t/h 2 )) with a 
decay time T2=h 2 /£ 2 r c : the fluctuation is smoothed out 
by the fast modifications of the reservoir configuration, 
and the line-shape becomes Lorentzian with a FWHM 
given by 2E 2 r c /fi.. Indeed, when motional narrowing oc- 
curs, the linewidth is reduced by a factor of the order of 
Et c /S. 

A comparison between our experimental data and the 
general expression of C (t) given above allows us to char- 
acterize the transition from the exponential to the Gaus- 
sian relaxation observed in Fig. ^ For this purpose, we 
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convolute Eq. QJ with the system response function (dot- 
ted lines in Fig. ^) which is obtained under white light 
illumination. The calculated fits are displayed in solid 
line in Fig. ^ We obtain an excellent agreement with in- 
creasing values of T,t c /H of 0.6 in (a), 1.05 in (b), and 1.35 
in (c), where the correlation time t c has a constant value 
of 10 ps. The gradual increase of T,t c /H demonstrates 
the transition from the fast modulation limit to the slow 
modulation one when increasing the incident power. As 
a matter of fact, it shows that motional narrowing oc- 
curs for low excitation. At first glance, one may think 
that increasing the number of photo-created carriers in 
the structure should induce some additional motion in 
the environment, so that the quantum dot line should be 
motionally narrowed in the high excitation limit. It is in 
fact the exact opposite, and we explain below the origin 
of this unconventional phenomenology. 



The existence of a fluctuating environment around the 
quantum dot originates from the presence of impurities, 
defects in the barrier material or in the wetting layer. 
Furthermore, this latter heterostructure is far from being 
an ideal two-dimensional quantum well so that the un- 
avoidable interface roughness gives also rise to localized 
states in the quantum dot surroundings 0, 0, . The 
characteristic parameters E and r c are thus determined 
by the population fluctuation in the various traps located 
around the quantum dot. We present Monte Carlo sim- 
ulations of the carrier dynamics in the reservoir for the 
understanding of the motional narrowing effect in quan- 
tum dots. 

In our model, the fluctuating reservoir consists in N 
uncorrelated traps which individually induce a Stark shift 
A of the quantum dot line. This simple scheme has been 
chosen for the sake of clarity. More realistic geometries 
can be used but do not alter the conclusions. Each trap 
is loaded with a characteristic capture time t\ , and emp- 
tied with an escape time t\ (Fig. Ufa)). In our time- 
dependent simulations, the random modifications of the 
traps population result in fluctuation over time SE(t) of 
the transition energy around a mean value Eq. Then, 
the spectral modulation amplitude E and the correlation 
time r c are calculated by computing the correlation func- 
tion of the energy fluctuation X(t), which is defined by 
the configuration average X(t) = (6E{t)6E(0)) 0. This 
correlation function is expected to have an exponential 
decay given by X{t)=Y? cy^{- t/r c ) 0. In Fig.^b), we 
display the correlation function X(t) computed for A^=50 
and averaged over one thousand configurations. We ob- 
serve on this semi-logarithmic plot the exponential decay 
of X(t), from which we extract E and r c . In Fig.^c) we 
display the correlation time r c normalized to the capture 
time Tj as a function of the ratio tj. /r-f . We observe that 
our calculations are fitted by using the simple expression: 



which shows that the correlation time is inversely pro- 
portional to the total fluctuation rate of the traps popu- 
lation. 

Finally, in Fig. EJd), we display the spectral modula- 
tion amplitude E normalized to its saturation value E s . 
Our numerical simulations are well reproduced by using 
the expression: 
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which can also be analytically derived by assuming ergod- 
icity of the reservoir. The saturation value E s of the spec- 
tral modulation amplitude is thus given by \/~NA/2. In 
semiconductor nanostructures, at low temperature car- 
rier escape is known to be much less efficient than the 
capture process 01- When tJt^<^1, the correlation 
time r c has approximately a constant value given by tj 
whereas E scales as 2E sAv /t|/t|. Therefore, an activation 
of the escape process by raising the incident power will 
lead to an increase of Er c /7i, in qualitative agreement 
with the phenomenology discussed above. 

On the basis of our simple microscopic model, we 
present a quantitative interpretation of our experiments 
at 10K as a function of incident power. In Fig. |3|we dis- 
play the measured values of the FWHM (squares) and 
Et c /S (circles, inset) as a function of the excitation den- 
sity, on a semi-logarithmic scale. These two parameters 
characterize the width of the line-profile and its shape 
since the crossover from the exponential to the Gaussian 
decoherence dynamics occurs around Er c /ft~l. We con- 
front our data with the theoretical values (solid lines) of 
the FWHM and Et c //i that are calculated with Eq. (|TI3j) . 
and we observe a fair agreement for the set of parame- 
ters E s — 400 ^eV, tj,~10 ps, and 1/t t = (1/t )\/P where 
tq~1.6 ns and P is in unit of 1 kW.cm~ 2 . 



The asymmetry between the power dependences of 
and Tf stems from the existence of different microscopic 
processes of various efficiencies in semiconductor physics 
[l2j. The constant value of 10 ps for Tj is consistent with 
an optical-phonon assisted capture in the traps around 
the quantum dot |15| . On the other hand, the inverse 
process of carrier escape by optical-phonon absorption 
is strongly inhibited at low temperature, so that it does 
not significantly contribute to the value of Tp. In fact, 
the escape rate dependence on incident power is char- 
acteristic of Auger-type processes. In the elastic colli- 
sion of two carriers where one is ejected from the trap 
while a delocalized carrier relaxes in energy, the escape 
rate is proportional to the density of delocalized carri- 
ers, which increases with incident power. In the expres- 
sion 1/t| = (1/t )v A P, the sublinear increase of the es- 
cape rate with P can come from the fact that the actual 
density of delocalized carriers is governed by bimolecular 
interband radiative recombination |l6j| or Auger scatter- 
ing |17| . which both result in a square root dependence 
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of the carrier density on P. Lastly, we note that in the 
limit of low incident power, the saturation value S s given 
by v^V A /2 may be power-dependent due to the reduced 
number of active traps that are populated from the reser- 
voir of delocalized carriers. Systematic measurements as 
a function of temperature and incident power on several 
quantum dots (to be described elsewhere) allow us to 
discriminate between the variations of N and tj in our 
experiments. The power-dependence of N is found to be 
negligible in the investigated range of incident power, in 
agreement with the assumptions made for our calcula- 
tions. 



Temperature-dependent experiments performed at low 
incident power (0.07 kW.cm~ 2 ) reveal the same striking 
phenomenology, namely we observe a gradual increase of 
Y>T c /h with temperature which demonstrates that mo- 
tional narrowing occurs at low temperature. However, 
the analysis of the data recorded as a function of tem- 
perature requires a more careful procedure because the 
line-profile exhibits broad side-bands around the central 
so-called zero-phonon line 0, . These side-bands are 
related to the radiative recombination assisted by the 
emission and the absorption of acoustic phonons. At low 
temperature (T<25K), their contribution to the total in- 
tegrated intensity of the photoluminescence signal is neg- 
ligible and can be safely neglected in the data interpreta- 
tion, as done for our first set of measurements as a func- 
tion of incident power at 10K. Above 25K, the presence 
of broad lateral side-bands around the zero-phonon line 
results in an interference contrast C (t) with an additional 
fast component, the decay of which follows our system re- 
sponse function (Fig. Ufa)). By adding a constant term 
to the line-profile extracted from Eq. QJ, we are able to 
fit the whole decay dynamics of the interference contrast 
(solid line, in green in (Fig. Ufa)), and thus separate the 
relative contributions of the lateral side-bands and zero- 
phonon line (line+empty circles, in blue in (Fig. ^a)). 

In Fig. ^Jb) we display the measured values of the 
FWHM (squares) and Sr c /fi (circles, inset) of the zero- 
phonon line as a function of the temperature. In 
the temperature-dependent measurements, the crossover 
from a Lorentzian to a Gaussian shape of the zero-phonon 
line occurs when increasing the temperature. More pre- 
cisely, a smooth transition analog to the one presented in 
Fig- fflis recorded in the temperature range between 30K 
and 45K. We then confront our data with the theoreti- 
cal values (solid lines) of the FWHM and T,r c /h that are 
calculated with Eq. (|1I3[) . by assuming that: 

- = -(l + n l (T)) + -(l + n 2 (T)) (4) 

Tl T\ T 2 

1 = L ni( T) + -n 2 (T) + - (5) 

Tf Tl T 2 T 3 

In Eq. © the constant term inversely proportional to 
T3 accounts for the contribution of Auger processes to 



carrier escape, as discussed above, and its magnitude is 
fixed by the incident power. In the interpretation of our 
temperature-dependent measurements, the relevant mi- 
croscopic processes for the carrier dynamics in the traps 
are phonon-assisted mechanisms where the capture (es- 
cape) is thermally-activated due to the emission (absorp- 
tion) of an acoustic phonon of mean energy E\ or an 
optical phonon of energy E 2 ■ The capture and escape ef- 
ficiencies are proportional to (l+ni(tj) and rn{t) respec- 
tively, where rii(t) is a Bose-Einstein occupation factor 
given by l/(exp(i^/£:T)-l). In Fig. 0{b), we observe a 
fair agreement between the data and the calculations for 
the set of parameters E s ~400 /ieV, ti~35 ns, T2~10 ps, 
r 3 ~6.5 ns, Ei^l meV, and E 2 ^20 meV. This latter ac- 
tivation energy is smaller than the characteristic optical 
phonon energy of 36 meV in bulk GaAs, and this differ- 
ence may come from strain and confinements effects in 
the traps formed in the GaAs barrier or in the InAs wet- 
ting layer |20j . As far as the carrier dynamics in the traps 
is concerned, we find that in our experiments performed 
below 50K, the capture process is dominated by optical 
phonon emission (t^t 2 ). On the other hand, the ther- 
mal activation of the escape one is predominantly due 
to acoustic phonon absorption below 40K, with an on- 
set of optical phonon absorption around 40K. Finally, we 
conclude that in the investigated range of temperature, 
we are always in the regime where the capture process is 
much more efficient than the escape one (t^/t-\<^1), in 
exact similarity to the power-dependent experiments at 
10K. 

The asymmetry of the capture and escape mechanisms 
is in fact the fundamental reason why motional narrow- 
ing strikingly occurs when decreasing the incident power 
or the temperature. If both processes had the same effi- 
ciency (TfWjJ, we would have r c ~T|/2 and which 
means that the spectral modulation amplitude would not 
depend on the time constant r^. Therefore, the ratio 
TiT c /h could only decrease when increasing the reservoir 
excitation. This situation corresponds to the well-known 
phenomenology from nuclear magnetic resonance where 
the activation of the nuclei motion induces the motional 
narrowing effect. In the present case where T|/t|<10 -2 , 
we are in the opposite regime where the correlation time 
is merely constant with relative variation smaller than 
10 -2 whereas the spectral modulation amplitude shows a 
steep increase with tJt^ (EcxyVjyTf ). The ratio Et c /S 
thus increases when increasing the reservoir excitation. 

We conclude that motional narrowing in the optical 
spectrum of a quantum dot broadened by spectral diffu- 
sion occurs in the unexpected regime of low temperature 
and low incident power. This result is of great impor- 
tance for the applications of quantum dots to quantum 
information processing. In fact, our study demonstrates 
that the decoherence dynamics can still be characterized 
by a decoherence time T 2 even in the presence of fluc- 
tuating local electric fields causing spectral diffusion. In 
that case, the zero-phonon line has a Lorentzian profile 
and its width is not given by the intrinsic radiative limit, 



■5 



but by the extrinsic reservoir fluctuation dynamics. In 
particular, we provide a straightforward interpretation to 
the widely debated issue of the origin of the temperature 
dependence of the zero-phonon linewidth. Eventually, 
we show that the motional narrowing effect strongly re- 
duces the environment-induced dephasing since the zero- 
phonon linewidth is approximately Sr c /ft smaller than 



the spectral modulation amplitude 1\j2ln1 S due to spec- 
tral diffusion. This unconventional motional narrowing 
opens a novel route towards the control of environment- 
induced dephasing for the applications of quantum dots 
to quantum information. 

*Electronic address: Guillaume.Cassabois@lpa.ens.fr 
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FIG. 2: (a) Schematic of the quantum dot and its fluctuating 
environment: rj. and t\ are the capture and escape times of 
carriers, (b) Correlation function of the energy fluctuation as 
a function of time, on a semi- logarithmic plot: Monte-Carlo 
simulations (squares) and exponential fit (solid line), (c) Cor- 
relation time r c normalized to the capture time t± versus rj,/r-f 
and (d) spectral modulation amplitude E normalized to its 
saturation value E c versus 7"j./rf: Monte-Carlo simulations 
(squares), and fits (solid lines) according to Eq. (|5J for (c) 
and Eq. © for (d). 




FIG. 3: Full width at half maximum (FWHM) and Er c /ft 
(inset) of the zero-phonon line of a single quantum dot, at 
10K. Data (symbols), and calculations (solid lines) based on 
Eq. Il l lot are plotted as a function of incident power, on a 
semi-logarithmic scale. 
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